© Ordnance Survey & Crown copyright and database right 2017 [add in AQMA shp file copyright notice]
© Ordnance Survey & Crown copyright and database right 2017 [add in AURN points copyright notice]
if(!is.null(nox)){
#plot(nox$date, nox$nox)
nox$day <- weekdays(as.Date(nox$date))
nox$month <- months(as.Date(nox$date))
nox$month<-as.ordered(nox$month)
ggplot(data = nox, aes(x = day, y=nox)) + geom_boxplot()
ggplot(data = nox, aes(x = month, y=nox)) + geom_boxplot()
} else {
print("No NOx measures")
}
## Warning: Removed 312 rows containing non-finite values (stat_boxplot).
# df <- get1Hdata('MY1', years=2015)
#
# # Aggregate to daily means and plot
# library('zoo')
# par(mai = c(0.5, 1, 0, 0))
# my1 <- zoo(x = df$Ozone, order.by = as.POSIXlt(df$datetime))
# plot(aggregate(my1, as.Date(as.POSIXlt(df$datetime)), mean),
# main = '', xlab = '', ylab = expression(paste('Ozone concentration [',
# mu, 'g/', m^3, ']')))